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Abstract. When subject to an external time periodic perturbation of frequency 
/, a Josephson-coupled two-state Bose-Einstein condensate responds with a 
constant chemical potential difference A/i = khf, where h is Planck's constant 
and k is an integer. We propose an experimental procedure to produce ac-driven 
atomic Josephson devices that may be used to define a standard of chemical 
potential. We investigate how to circumvent some of the specific problems derived 
from the present lack of advanced atom circuit technology. We include the effect 
of dissipation due to quasiparticles, which is essential to help the system relax 
towards the exact Shapiro resonance, and set limits to the range of values which 
the various physical quantities must have in order to achieve a stable and accurate 
chemical potential difference between the macroscopic condensates. 



PACS numbers: 03.75.Lm, 74.50.+r, 05.30.Jp, 06.20.Fn. 

1. Introduction 

The realization of Bose-Einstein condensation in dilute atomic gases jll2!3j has opened 
the possibility of investigating the macroscopic behaviour of coherent quantum matter 
in a new class of physical systems different from superconductors and helium liquids. 
Together with vortex formation, the Josephson effect between weakly coupled Bosc 
condensates is probably the most characteristic signature of superfluidity. Initially 
predicted 0] and observed [HUE] in the context of superconductivity the Josephson 
effect has also been observed in superfluid 3 He 0|Ej, a fact that underlies the 
fundamental nature of quantum behaviour at the macroscopic scale. Atomic Bose- 
Einstein condensates (BECs) have it in common with superconductors and superfluids 
that their most remarkable properties are derived from the existence of a macroscopic 
wave function. Because of the profound analogies between these gauge-symmetry 
broken systems, the possibility of observing the Josephson effect in trapped atomic 
gases was early recognized [§l ll()llT11 li2|. Although some preliminary evidence for the 
existence of the Josephson effect in atomic gases already exists |13U14) , it is fair to say 
that, compared with its superconducting and superfluid counterparts, the exploration 
of the physics of weakly linked atomic condensates is still in its infancy. Given the 
potential richness and the convenient tunability of BEC setups, it is clear that decided 
progress must be made to investigate this whole new class of problems. At present 
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however, an important limitation is that the technology of atom circuits (also called 
"atomtronics" ) is still moderately developed, although a bright future can already be 
foreseen |15I16| . In this sense, proposals that do not rely on the feasibility of controlled 
atom transport are particularly welcome. 

One of the most remarkable properties of superconducting devices is the existence 
of Shapiro resonances. When two weakly coupled superconductors are subject to 
a voltage difference that is the sum of a dc component V and a periodic signal 
u sin(2-7r/£), a continuous range of nonzero dc currents are possible if 



where 2e is the Cooper pair charge, h is Planck's constant, and k is an integer |rjlll7|. 
A more common setup is one in which, due to the existence of a large impedance 
in series, an external current source is applied. If the external current consists 
of a dc contribution and a weak ac perturbation of frequency /, the Josephson 
link displays a dc V-I characteristic with voltage plateaus at values satisfying the 
resonance condition QJ. The height of the Shapiro steps linking the voltage plateaus 
provides a method to measure the constant of nature 2e/h with such precision and 
universality (181119) that, since 1972, the reversed view has been adopted whereby 
2e/h is assumed to be known and Equation is used to define a standard unit of 

voltage [niEniEi]- 

In this paper we address the question of whether it is possible to prepare a BEC 
Josephson junction (BJJ) that satisfies the resonance condition JTJ in a stable form. 
Another way of phrasing the problem is whether macroscopic self-trapping of a large 
population imbalance 22, 22] may be stabilized against dissipation by an oscillating 
perturbation. A condensate in a time-dependent trap |24j . and particularly a double 
condensate [21], has been shown to provide a convenient test ground for quantum 
chaos. Here we face a different regime which predominantly involves regular motion 
with dissipation. The challenge is to identify the conditions under which dissipation 
can be exploited to let the system evolve towards one of the stable resonance islands. 

At present we have a limited understanding of the dissipation mechanisms 
operating in a BJJ which render its macroscopic phase-number dynamics non- 
conservative. We know that contributions may arise at least from incoherent exchange 
of normal atoms |llj , creation of quasiparticles by the fluctuating condensate |23j , or 
spontaneous atom losses |26I27) . As long as the net atom loss stays relatively small, its 
effect on the phase dynamics is similar in many respects to that of incoherent particle 
exchange, both yielding Ohmic dissipation under a wide range of circumstances. 

We wish to stress that the physics discussed in this paper applies both to double- 
well condensates, displaying the external Josephson effect, and to optically coupled 
two-component atomic condensates, which exhibit the internal Josephson effect |12| . 
A central concept in the forthcoming discussion is the application of a time-dependent 
external potential. This may be achieved by the application of suitably designed 
time-dependent magnetic fields or dipole forces. Another essential ingredient, not 
considered in Ref. |22]> is that the Josephson current can counteract the dissipative 
current, thereby permitting the existence of stationary Shapiro resonances. Moreover, 
by providing an effective friction to the phase dynamics, the dissipative current ensures 
that the system gets arbitrarily close to the exact resonance at sufficiently long times. 

We want to remark that no fundamental reason prevents the possibility of 
exploring this novel approach to Shapiro resonances in the context of superconductors. 
There, capacitive couplings might be used to control or measure voltage differences. 



V = 




(1) 
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Section |21 is devoted to a quantitative formulation of the problem. We know that 
each phase-locked solution generates an attractive basin in phase space. We raise the 
question of how dissipation can be exploited to let the system land successfully on 
that region never to escape from it except for quantum tunnelling decay. The solution 
is offered in Section where we investigate the dissipative dynamics in phase space. 
Later, in Section 01 we provide an estimate of the probability that a system prepared 
with a large chemical potential difference evolves in such a way that, in the process 
of decaying, it is trapped in a classically stable well. The stochasticity comes from 
the intrinsically random nature of the choice of phase made by the double condensate 
system shortly after the connection [221 ■ Section |5l is devoted to the effect of quantum 
fluctuations, which may cause macroscopic tunnel decay. The combined analysis of 
Sections EHHl yields a set of constraints that must be fulfilled to achieve an optimal 
realization of the predicted Shapiro resonances. A concluding discussion is provided 
in Sectional 



2. Formulation of the problem 

We consider N Bose-condensed atoms in a double- well trapping potential V(x, £) that 
depends on a harmonically time-dependent control parameter = £ocos(fit). In a 
situation sufficiently close to equilibrium, a split condensate is formed which can be 
described within a two-mode approximation. Thus, we use the for the field operator 
the ansatz if>{x) = ifA(x)a + (fB(x)b, where the orbital functions (fA,B are normalized 
to unity. This yields the two-mode Hamiltonian |1(JII12| 

H = 6 + tfa) + E A (N A ,{) + E B {N B , £)• (2) 

Here, lor denotes the effective Rabi frequency of the two-mode problem, which is 
proportional to the Josephson coupling energy |12j . It is important to note that, in 
the case of optically coupled two-component BECs (internal Josephson effect), lur 
does not generally coincide with the Rabi frequency fl r governing atomic transitions. 
Rather, one has ujr = £IrSab, where sab = Jdx(p A ipB EH- It is only when the 
trapping configuration and the interactions are such that (Pa(x) = </?s(x), that the 
two frequencies become identical. The Gross-Pitaevskii energy of fragment i, 

EiiN^O = Ni Jdx<p*(x) (-¥-A + V(x,0 + ^|^(*)| 2 ) <Pi(x), (3) 

i = A, B, has inherited a time-dependence from the control-parameter 

For an analysis of the system in its classical limit, we replace in the Heisenberg 
equations of motion the operators a, b by \ / NAe~ t<f ' A , ^fNBe~ l<t ' B and expand the 
Gross-Pitaevskii energies for small £ and small number imbalance compared with 
the equilibrium values N A and Ng (usually taken to be N/2). For the variables 
<fi = (<pA — 4*b) and z = (Na — Nb)/N = 2n/N, with time measured in units of 1/lur, 
this yields 

Z 

4> = —j=^ cos0 + Az + ecos(^) = Afi/h, (4) 

z= — V 1 — z 2 sin tj) — 70, (5) 
with the scaled driving amplitude 



go 



d 2 E A (N A ,Q d 2 E B (N B ,Q 
d£,dN A d£,dN B 



(6) 

N A =N B = N/2,£=0 
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Figure 1. Phase space portrait of the undriven [e = (a)] and the driven [e = 100 
(b)] undamped momentum-shortened pendulum Hamiltonian JSJ at stroboscopic 
times tmod27rAi/f2 = tt/2. The interaction is A = 10 4 and the driving frequency 

n = iooo. 



and the effective interaction constant 



A 



N 
2hw, 



d 2 E A (N A ,0) d 2 E B {N B ,0) 



dN% 



(7) 



N A =N B =N/2 

Equations Q, J5J can also be obtained from the classical non-rigid pendulum 
Hamiltonian 10, 12 (hereafter energies are expressed in units of Itwr) 

H(z,(j),X) = - y/l - z 2 cos <f) + ^Az 2 + ezcos(ilt), (8) 

where (z, (f>) are canonically conjugate coordinates. For a simplification, we have 
assumed a symmetric situation with Ea(N/2, £) = Eb(N/2, — £) — the generalization 
is of course straightforward and results in an additional phase drift. The last term 
in the equation of motion has been introduced phenomenologically to describe a 
dissipative current, i.e. an incoherent exchange of atoms. For high temperatures, 
ksT 3> A/i, this current is Ohmic, i.e. proportional to the chemical potential 
difference, h = -GA/i 11 . Following the reasoning by Josephson 01, the chemical 
potential difference is given by the time derivative of the relative phase, </> = A/i/h, 
and thus we obtain the dissipative term in equation (J3J with 7 = 2HG/N |. We will 
also analyze the more general case in which n may not be a linear function of <j), a 
situation which is likely to appear at low temperatures. 

In order to put the central discussion into the proper context, it is convenient 
to review the dynamics of the relative phase after two independently prepared 
condensates are connected [22] • Before the connection, the phase is completely 
undefined, which means that the Fock state of fixed particle number is in a coherent 
superposition of different phase states. Upon connection, various mechanisms 
involving quasiparticle dynamics intervene to destroy the coherence between the 
different phase states. As a result, the reduced density matrix of macroscopic phase- 
number system becomes quasi-diagonal in the phase representation, which is to say 

I At first sight there seems to be an ambiguity as to whether dissipation should come as — 7</> or, 
rather, —-y'z. That the former is the correct approach can also be inferred from a careful study of 
quantum dissipation models 1291 . In the present context, however, the choice is without practical 
consequences. 
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Figure 2. Population imbalance z(t) for interaction A = 10 4 , damping 
7 = 3 X 10 , driving amplitude e = 100, and frequency Q = 1000 for two different 
initial values of the phase. The saturation time-averaged value z = 0.1 of the 
upper solution corresponds to the chemical potential difference Ml (k = 1). The 
inset is a blow-up that resolves the fast underlying Bloch (or MQST) oscillations. 

that the phase is effectively measured among a menu of values uniformly distributed 
between and 2tt. After this quick definition of the phase, a phase-number Gaussian 
wave packet forms that evolves semiclassically in the parabolic tight-binding lattice 
formed by the different number eigenstates. Due to the interaction term, the system 
experiences a tilted lattice and, thus, undergoes Bloch oscillations, which may be 
viewed as the ac Josephson response to an approximately constant A/x created by 
the number imbalance with the possible concurrence of other factors. Its trajectory 
follows one of the running solutions in phase space diagram shown in Figure^. Thus, 
the system oscillates around a nonzero number average z displaying the so-called 
macroscopic quantum self-trapping (MQST) However, these oscillations are not 
stable, since the fluctuating condensate atom number in each well couples to the 
many-quasiparticle field and experiences dissipation. Thus the energy stored in the 
macroscopic degree of freedom decreases and z decays slowly towards its equilibrium 
value [2"3| . 

In the presence of an ac driving, the situation may change qualitatively, since 
classically stable resonance islands form in phase space for values of z, <f> such that 

Az = kto, dz/dt = 0, (9) 

where k is an integer (see Figure H 3 ). Here, the bars indicate the time average over 
fast Bloch or external driving oscillations, of frequencies fcO and f2, respectively, and 
the time derivative is meant over a longer time scale. In general, the system displays 
running solutions weakly oscillating around an average value z which slowly decays 
because of dissipation. Such a decaying trajectory may or may not be trapped in one 
of the attractive basins around the Shapiro resonances. 

For a provisional answer to the central question formulated above, we integrate 
numerically the equations of motion and JSJ in the Josephson regime (A ^> 1; 
we take A = 10 4 ), where the BEC Josephson junction displays, in the absence of 
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dissipation, stable running MQST solutions. When starting with an initial imbalance 
z(0) that lies above a resonance island, we find two qualitatively different types of 
solutions: The first type (lower curve in Figure^ is what one would typically expect, 
namely that, owing to the Ohmic current, the imbalance decays until it leaves the 
MQST regime and ultimately exhibits damped plasma oscillations. Except for the 
small kink displayed when crossing the resonance zone, the dynamics is qualitatively 
the same as in the undriven case. However, for the same parameters but a different 
initial phase, there exist also solutions with an intriguing behaviour (see upper curve 
in Figure^: There, the system gets trapped in the resonance island and the transient 
decay comes to a standstill instead of continuing towards the true equilibrium solution 
[z = 0). At long times, the number imbalance z oscillates with amplitude e around 
a non-zero value z\ = fl/A. We will show that then the average chemical potential 
difference (in scaled units) is precisely a multiple of the driving frequency. More 
specifically, we find 

Afi = kQ + ecos(nt). (10) 

We wish to remark that the BEC setup considered here differs substantially form 
its superconducting analogs in that there the stabilization of a voltage difference at 
the Shapiro resonance is achieved with the concourse of an externally imposed current, 
i.e., in an experiment a current source is needed. Here by contrast, the stabilization of 
a chemical-potential difference cannot rely on such external current sources which for 
BECs are not (yet) available. As a drawback, the present type of "voltage" standard 
can be obtained only at the price that a single run of the experiment does not warrant a 
stable solution. In the following, we provide a detailed analysis of this type of Shapiro 
resonances, estimating their parameter dependence and, eventually, the probability 
that they can be spontaneously realized. 



3. Realization of phase-locked solutions 

Our central interest lies in the existence and the stability of so-called phase-locked 
solutions, i.e., solutions (j>{t) that have only small fluctuations around a long-time 
average (f> w kflt. Then, since the main chemical-potential difference between the 
two condensate fragments comes from the interaction energy, the average number 
imbalance z settles down at a value Zk = k£l/A. This motivates the ansatz 

= kQt + ^ sin(fit) + 84>, (11) 
fcf2 

z = — + a cos(fcfii) + Sz, (12) 

where a denotes the yet unknown amplitude of the fast residual Bloch oscillations 
that can be appreciated in the inset of Figure El We would like to choose a such that 
the resulting 5(f) and Sz display only slow, decaying oscillations around a stationary 
solution 64> = const., Sz = to which they tend at long times. In other words, a should 
be such that, in the absence of Sz and with 8cj> replaced by a constant, equations 
<(T2"|> correctly capture the asymptotic system dynamics. 

The equations of motion for 8<p and Sz are obtained by inserting Q into 

the original equations JHJ. The resulting system of equations involves rapidly 
oscillating coefficients of period 27r/f2. It can be shown that, if a = l/k&l and fl 1, 
then 5<j) and Sz vary more slowly than those coefficients (a result which we have 
confirmed by numerical studies), so that the different time-scales can be separated 
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Figure 3. Sketch of the effective washboard potential 1151 as a function of the 
phase difference Srp for ■ykQ < Jf c (e/Q). 



consistently and all time-dependent coefficients can be replaced by their time-averages. 
On the other hand, we are interested in the case of two condensate fragments which 
are comparable in size (\z\ -C f). Thus we neglect within this analytic discussion 
the momentum shortening and replace the square root in the Hamiltonian JSJ) by 
unity |12| . Finally we obtain the equations of motion 

±6cf> = A6z, (13) 

^_S<f> = - Jk(e/fl) sm(Sct>) - jktt - jASz, (14) 

which describe a dissipative particle in the static tilted washboard potential sketched 
in Figure They can be obtained from the Hamiltonian 

H(5<p,6z) = ^ASz 2 - J k (e/n)cos(S<p) +-fknS(P (15) 

together with the dissipative force Fdiss = —jASz. The new canonical coordinates Sz 
and 5(f> represent momentum and position, respectively, and Jk denotes the fc-th Bessel 
function of the first kind. We emphasize two differences between these equations and 
the undriven rigid pendulum Hamiltonian: First, there is a tilt ~/k£l which originates 
from a constant dissipative current caused by a finite number imbalance. Despite its 
physical origin, it appears formally as a conservative force. Second, the Josephson 
coupling energy for the phase shift 8<j) is renormalized by a factor Jk(s/fl) which, for 
e <C f^, renders it much smaller than the original Josephson energy. Correspondingly, 
the plasma frequency is renormalized by a factor *J JfJe/TT)§. Once S(p is trapped 
within one specific well, it exhibits damped plasma oscillations with the renormalized 
frequency until it ultimately comes to rest. In this stationary solution one has 
Scf> — Sz — 0, which is possible thanks to the cancellation of the first two terms in the 
r.h.s. of Equation 1)14(1. Therefore, one finds from Equation 1)11(1 and the Josephson 
relation <j> — A/J./H the chemical potential difference ((1011 . 

§ If we had relaxed the assumption \z\ 1, then Jj.(e/£7) would have appeared multiplied by a 
factor yi — z 2 . However, one must remember that, if |z| were to become comparable to unity, then 
the two-mode approximation might have to be revised. 
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Having mapped the long-time dynamics of the originally time-dependent situation 
to an equivalent static problem, we are now able to derive the conditions under 
which phase- locked solutions exist. This results in the restrictions of the parameters 
summarized in Table ^ Obviously, the effective washboard potential in the 
Hamiltonian i|15|) possesses minima only for jkil < Jfc(e/Sl) which amounts for small 
driving amplitudes e <C SI to -fktt <C (e/2il) k . This is most easily satisfied for k = 1, 
where the condition for the existence of wells reads e > 27SI 2 . Since on the other 
hand, the driving should be weak, e <C SI, both conditions on the driving amplitude e 
can be fulfilled simultaneously only for sufficiently small dissipation, 7 <C I/O. 

Also the interaction strength A obeys restrictions, which are obtained in the 
following way: Physically, the idea behind our scheme is to counterbalance a dissipative 
current with the help of ac driving, thereby stabilizing a MQST solution against 
friction. The condition for operating clearly in a regime where MQST dominates is 
2 / vA~ <c z |22I23| . On the other hand, the two condensate fragments should not differ 
too much in size, thus, z <C 1. Since we aim at stabilizing the imbalance z = SI/ A, 
fulfilling both conditions requires a sufficiently large interaction strength, typically 
A !()''. 

A further physical reason for operating the Bose-Josephson junction in the 
(interaction-dominated) Josephson limit stems from an important property of phase- 
locked solutions, namely, that generally the centre of an attractive basin is not 
truly at rest, but follows the trajectory of a particle in the absence of driving and 
dissipation [SOI - In the Josephson regime, such a behaviour corresponds to the MQST 
solutions exhibiting Bloch oscillations with amplitude a — 1/Az around an average 
value z, as can be appreciated in Equation H12(l . At resonance, their relative amplitude 
becomes a/z = 1/Az 2 = A/57 2 . Ideally, for the realization of stable phase- locked 
solutions, this amplitude should be small, as is the case in the Josephson regime 
(A ^> f). On the contrary in the (non- interacting) Rabi regime (A < 1) 12 , the 
undamped and undriven pendulum exhibits large number (Rabi) oscillations around 
z = 0, i.e. there is no MQST. This results in large oscillations of any attractive basin 
which render BJJs in the Rabi regime useless for the present purposes. 

We conclude this section with a comparison between Shapiro resonances and 
another well-known, apparently similar dynamical phenomenon. It has been shown 
[21} that for a high-frequency driving with SI ^> A^i/h, the Josephson energy acquires 
a factor Jo(e/£l) and that, consequently, the coherent exchange of atoms is brought 
to a standstill if e/fi is chosen as a zero of the Bessel function Jo- Such an effect 
amounts to an interacting version of the so-called coherent destruction of tunnelling 
(CDT) predicted for single particles in bistable potentials [32 . We emphasize that the 
present phenomenon is different in two respects: First, CDT takes place at zeros of 
Bessel functions and, thus, requires large driving amplitudes, while here, the driving 
amplitude is much lower. Second, CDT fades out under the influence of dissipation, 
while in our case, moderate dissipation is essential for the convergence towards a 
phase-locked solution. 

4. Locking probability 

After investigating which is the parameter regime that permits stable phase-locked 
solutions, we turn to a remaining intriguing question: What is the probability that 
such a solution is hit by starting with a random initial phase? To give a crude estimate 
for the answer, we employ again the analogy to the dissipative motion in a tilted 
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Table 1. Conditions for the existence of phase-locked solutions which are suitable 
for a chemical-potential standard with Afi = Ml (i.e. z = Q/A and k = 1). 

restriction physical significance 

e -C Q chemical potential difference much larger than its modulation 

e > existence of stable wells in the washboard potential 

A <g f2 2 operation within MQST regime, Bloch oscillations small 

A>!1 number imbalance small 



washboard potential. 

Let us assume that we start out of resonance with an imbalance z\ < z < Z2, 
with Zfe = &S7/A, that corresponds to Sz > 0, i.e., to a particle moving uphill in the 
tilted washboard of Figure |3| Having only finite energy, the particle will bounce in 
one specific well. At the entry point, the energy will be in the range [E^, E^ + 27T7fcS7], 
where E\, is the potential energy at the top of the lower barrier. While moving within 
the well during one cycle, the particle dissipates the energy 

£diss W - / ^diss # « 27T 7 AV2£/AO. (f 6) 

Here, we have estimated the maximum value of the dissipative force -Fdiss = — 7A Sz 
[cf. Equation l|14|) ] from the maximum kinetic energy ^A5z 2 w 2Jfc(e/S7) ss e/ft for 
k = 1 and 7 <C e/Sl 2 . Assuming for the particle's phase space trajectory the shape of 
an ellipse, yields the r.h.s. of Equation If the initial energy minus the dissipated 
energy lies below the barrier, the particle will end up at rest in the well. By assuming 
that the random initial phases translate into equally distributed initial energies, we 
find that this happens with probability 

^^"V"^' (17) 

if -Ediss < 27T7f2, and w = 1 otherwise. Ideally, one would like to have 10 = 1 to ensure 
phase locking in all runs. However, it is easy to show from the constraints in Tabled 
that this would require fi<£, which, for reasons already indicated, is not of physical 
interest. Thus one must content oneself with identifying a range of parameters that 
make the probability of landing in a resonance non-negligible. 

One might conceive situations where the control of the relative particle number 
before the connection (as well as of the interaction parameters) were so good that 
the system could be forced to be within the attractive basin from the start. While 
designing z(0) to be very close to the resonance value z = 11/ A would surely increase 
the probabilities of relaxing towards the Shapiro resonance, we wish to emphasize here 
that success can never be fully guaranteed because of the intrinsically random nature 
of the initial phase <fi(0), as can be clearly inferred from inspection of Figure 

To confirm the analytical estimate (| 1 7|1 for the locking probability, we have 
integrated numerically the equations of motion (0J and JSJ starting with a number 
imbalance z(0) = 1.5J7/A, and a random initial phase </>(0). The parameters have 
always been chosen such that e <§; SI A and 7 <C according to the requirements 
derived in the last section. Figure 0] compares the analytical estimate (|17fl to the 
fraction of numerical runs that converge to a phase-locked solution. The numerical 
and the analytical result agree well in the regime SI > 30. For SI = 100, the locking 
probability typically assumes values of order 10%. We attribute the worse agreement 
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Figure 4. Estimate 1171 for the locking probability compared to the fraction 
of phase-locked solutions (symbols) obtained from numerical integrations of Eqs. 
Ill and J5J with random initial phase. The driving amplitude is e = 0.1 f2 (full 
line, filled symbols) and e = 0.01 Q (broken line, open symbols). The interaction 
is A = Q/z-l where z\ = 0.1, and the dissipation 7 = 0.1 e/f2 2 (circles) and 
7 = 0.05 e/C 2 (triangles). Each data point is obtained from 1000 simulation runs. 



between theory and simulation for Q < 30 to the fact that, in such a range, fl becomes 
so small that the requirement A <C fi 2 (see Table can no longer be satisfied if, as is 
the case in Figure 01 A is constrained to be fl/zi with z\ = 0.1. An extreme case of 
disagreement is found for the lowest frequency considered when e = 0.1 f2. 

The Ohmic dissipation which we have assumed in our previous analysis is only 
justified for temperatures well above the chemical potential difference |11|. For 
a typical condensate however, the chemical potential is of the same order as the 
temperature and, thus, the dissipative current may not follow a linear law. Therefore, 
as a last item within the classical analysis, we consider the more general dissipative 
current i|diss = ~g( z ): which might depend on various parameters like, e.g., the 
interaction strength and the temperature. Using the ansatz (I12|) . we linearize 
g(z) w —g(zk) — g'(zk) 5z and repeat the analysis from above. We find for k — 1 again 
the locking probability 117|) but with an additional factor z\g' {z{) j g(z\) = z\ [In g(zi)}' , 
z\ = CI/ A. On the other hand, the slope of the washboard, i.e. the average dissipative 
current, acquires a factor g(z\). For very low temperatures ksT A/x, for instance, 
the quasiparticle-excitation decay mechanism discussed in Ref. |23j yields g(z) oc z 2 , 
and therefore the locking probability is enhanced by a factor of 2, but still is of the 
same order. 

5. Quantum fluctuations 

The classical treatment of the Bose-Josephson junction given above, has its limitations 
due to quantum fluctuations of the number n and the phase <j). They come from the 
commutation relation [n, e 1 ^] = e 1 ^ which amounts for small phase uncertainty to 




2i/N. 



(18) 
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This leads basically to two constraints for the washboard potential: First, the potential 
wells of a washboard potential are quantum mechanically metastable since a particle 
will tunnel out at a rate \6'6\ 

n = uj e- 2nEo/hul0 . (19) 

To justify the classical treatment from above, 1/k must be larger than all other 
time scales of the problem. A second, related point is that the phase-space region 
corresponding to the metastable well must support sufficiently many quantum states 
to treat both the number and the phase as continuous classical variables. For a well 
of depth E with curvature Wq, the number of quantum states can be estimated as 
r7i = Eq/Hluq. From both arguments, we conclude that the ratio of the potential depth 
and the energy quantum of a small oscillation at the bottom of the well determines 
whether we operate in the classical limit Eq ^> ItjjJq. 

To determine the connection between Eq, Hlvq and the non-standard commutation 
relation l|18l) and our scaled parameters for k = 1, we have to accomplish the 
replacements 

E ^-, h^-, wg-^. (20) 

This results in 



S^vSi' (21) 

which is essentially the ratio between the Josephson coupling energy and the 
renormalized plasma frequency in the metastable well. For typical parameters used 
above (e — 0.0151, A = 10 O, O > 100), a condensate consisting of N > 10 5 atoms 
supports m « 100 states and yields an escape rate that is practically zero. Therefore, 
we do not expect any relevant quantum correction. 



6. Conclusions 



We have investigated a realistic setup that may provide the basis for a standard 
of chemical potential difference between weakly connected atomic Bose-Einstein 
condensates. Due to the still rudimentary development of atom circuit technology, 
we have focused on schemes that do not require the coupling to a an external circuit, 
thus staying away from straightforward analogues of well-tested superconducting 
devices. In particular, we have investigated the possibility of connecting two separate 
condensates in the presence of ac driving in such a way that, as the system relaxes 
towards equilibrium, it has an appreciable probability of being trapped in a Shapiro 
resonance for which A/i is exactly an integer multiple of hf, with / the driving 
frequency. If the resolution of the imaging process is good enough, a frequency 
= 2itf can be found satisfying simultaneously 

AAz<0<A. (22) 

While the second inequality ensures a small number imbalance (see Table ^ , the 
first one expresses the experimental ability to resolve two different Shapiro plateaus, 
Az < Zk+i — Zk = f2/A. Therefore, under any reasonable visibility conditions, it 
should be possible to identify the integer ratio A/i//i/, i.e., we can know the precise 
Shapiro step where the double condensate has become locked. Thus, we see that the 
main requirement of a precision measurement is satisfied, namely, that, with the help 
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of an accurately controlled frequency, a poor measurement of a physical quantity (the 
atom number imbalance) enables a fine measurement of another physical quantity (the 
chemical potential difference). 

Instead of being a hindrance, dissipation provides here a crucial help, since 
it allows the system to spontaneously relax towards a stable resonance where the 
chemical potential difference is guaranteed to have a precise value. We have identified 
a number of parameter constraints that must be satisfied in order to have an optimum 
control of the Shapiro resonance. These include the need to have a weak driving signal, 
the operation within the collective Josephson regime, the need to keep average number 
imbalance small, the requirement of even smaller Bloch oscillations, the necessity to 
optimize the probability of landing in the desired Shapiro step out of a random initial 
phase, and the robustness against escape by quantum tunnelling from the metastable 
state. 

Fortunately, this set of constraints does still allow for a window of realistic 
parameters within which accurate chemical potential differences could be realized. 
Let us for instance consider N = 10 6 23 Na atoms in a split trap with frequency 
Who = 2-7T x 100 Hz. The bulk properties of the condensate are readily estimated within 
the Thomas- Fermi approximation to yield the chemical potential /i sa 70 thuj^o and the 
mean field energy Ea(N) HN 7 ^ 5 x 0.01 s _1 . Estimating the Josephson coupling 
energy, respectively the Rabi frequency, is more tedious and relies sensitively on the 
assumptions made for the shape and the size of the barrier. We assume here a Gaussian 
shape of width 6 firm and height of 1.05 /i such that for a driving amplitude e = 0.01 CI, 
the top of the barrier lies always above the instantaneous chemical potential of both 
wells. Along the lines of Refs. |11II12| . we obtain the effective action So = 1.7 h and, 
thus, the Rabi frequency u>r w 2ir x 0.05 Hz. This finally yields the effective interaction 
A w 10 4 that we have assumed above. The (bare) Josephson plasma frequency is 
approximately u>jp = 2ir x 5 Hz; for the resonant solution in Figure [3] phase locking 
sets in after 5 s. 

One might wonder what may be the use of a "voltage" standard in the absence 
of an atom circuit. First of all, no fundamental reason prevents the use of the present 
scheme in future atom circuits. Then, this concept will have to be tested against ideas 
more directly borrowed from the technology of superconductors. Second, already 
within the currently limited possibilities of condensate transport, the realization of a 
Shapiro resonance may provide a convenient playground for the investigation of novel 
BEC scenarios. For instance, it seems possible to extend the experimental work of 
Ref. m] to include an external ac driving that locks a pair of wells into resonance. 
The static behaviour of that pair of condensates should contrast with the dynamic 
behaviour of atoms in neighbouring wells. That would provide a most direct test 
of the Josephson effect in atomic condensates. Third, a fine control of the chemical 
potential may open the possibility of detailed checks of our current understanding 
of the novel many-body problem posed by trapped quantum gases. Quite likely, the 
effect of interactions will be easily isolated from that of gravity or other well-controlled 
external fields. Finally, there is the hope that the present work may provide the basis 
for more efficient or practical concepts that will make the emergent "atomtronics" a 
precise technology. 
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